Multilevel effects in the Rabi oscillations of a Josephson phase qubit 
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We present Rabi oscillation measurements of a Nb/A10 x /Nb dc superconducting quantum in- 
terference device (SQUID) phase qubit with a 100 /im 2 area junction acquired over a range of 
microwave drive power and frequency detuning. Given the slightly anharmonic level structure of 
the device, several excited states play an important role in the qubit dynamics, particularly at high 
power. To investigate the effects of these levels, multiphoton Rabi oscillations were monitored by 
measuring the tunneling escape rate of the device to the voltage state, which is particularly sensi- 
tive to excited state population. We compare the observed oscillation frequencies with a simplified 
model constructed from the full phase qubit Hamiltonian and also compare time-dependent escape 
rate measurements with a more complete density-matrix simulation. Good quantitative agreement 
is found between the data and simulations, allowing us to identify a shift in resonance (analogous to 
the ac Stark effect), a suppression of the Rabi frequency, and leakage to the higher excited states. 

PACS numbers: 74.50.+r, 03.67.Lx, 85.25.Dq 



I. INTRODUCTION 

Over the past decade, a variety of qubits based on 
superconducting Josephson junctions has been proposed 
and experimentally realized.— i 2 -^^ Among them is the 
phase qubit, which can take the form of a single current- 
biased junction, 6 '- a flux-biased rf superconducting quan- 
tum interference device (SQUID) j£ a low inductance dc 
SQUID j2 a large inductance dc SQUID where one of the 
junctions is thought of as the qubit fiS or a vortex in 
an annular junction^ The behaviors of the devices are 
similar, as they all employ a large area junction, where 
the dynamics of the quantum mechanical phase differ- 
ence across the junction (or the orientation of the vortex 
in the last case) are determined by a tilted washboard- 
like potential. The two lowest states in a well of this 
potential serve as the qubit basis. Among the possible 
benefits of the phase qubit are its relative insensitivity 
to charge and flux noise and its ability to operate over 
a wide range of parameters. Recent demonstrations in- 
clude fast readout^ Rabi oscillations^ ' 10 ' 12 ' 13 and sim- 
ple capacitive coupling of two qubits as measured through 
spectroscopy^ simultaneous state measurement^ and 
state tomography.— 

In the phase qubit, higher excited states can impact 
the dynamics strongly, due to the nearly harmonic level 
structure within the potential well. The presence of these 
quantized energy levels has been clearly detected by mon- 
itoring the decay of a thermal population^ 7 and by mi- 
crowave spectroscopy with single^ and multiphoto n 12 ' 19 
transitions. Multilevel systems can, for example, be ex- 
ploited for state initialization^ and readout^ quantum 
logic gate a 21 ' 22 ' 23 and algorithms^ and cryptography. 25 
However, when controlling the state of the phase qubit 
with a microwave current, off-resonant excitation of 
higher levels leads to leakage out of the desired qubit 



space and therefore a loss of quantum information. 



Rabi oscillations serve as both a standard demonstra- 
tion of quantum state manipulation and a diagnostic for 
decoherence and control fidelity. While these oscillations 
take on a simple form in a two-level system, several ef- 
fects can occur in a multilevel system: the oscillations 
can distort, the higher states can become populated, and 
the basic resonance properties of the oscillations can un- 
dergo subtle shifts. Ultimately, to achieve fast and accu- 
rate control of the qubit state, all of these effects need 
to be carefully characterized. Progress in this area has 
been predominantly theoretical in nature, pointing out 
the conditions under which errors are introduced and 
methods to minimize their effects ! 12 ' 26 ' 27 ' 28 ' 29 ' 30 There 
has, however, been less direct experimental evidence con- 
firming the validity of models upon which these predic- 
tions are based ) 9 ' 12 ' 31 ' 32 particularly with regard to leak- 
age. 



In this paper, we investigate the influence of the higher 
states of a phase qubit by examining the behavior of 
Rabi oscillations at high power, where their impact is 
greatest. In Sec. lU we describe the design of our qubit 
and a Hamiltonian that approximates its dynamics. Also 
discussed are various experimental details and the read- 
out scheme that allows measurement of very small ex- 
cited state populations. Section UlTl describes Rabi oscil- 
lations taken at a range of power and detuning from res- 
onance, along with comparisons to a simple theory, while 
in Sec. IIVI we develop a more complete density-matrix 
model including the effects of decoherence and noise. Fi- 
nally, Sec. fVl contains a summary of the key points of our 
findings. 
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FIG. 1: The dc SQUID phase qubit. (a) The qubit junction 
Jl (with critical current Ioi and capacitance Ci) is isolated 
from the current bias leads by an auxiliary junction J2 (with 
I02 and C2) and geometrical inductances Li and L2. The 
device is controlled with a current bias If, and a flux current 
If which generates flux <& a through mutual inductance M. 
Transitions can be induced by a microwave current I r f, which 
is coupled to Jl via C r f. (b) When biased appropriately, the 
dynamics of the phase difference 71 across the qubit junction 
are analogous to those of a ball in a one-dimensional tilted 
washboard potential U. The metastable state \n) differs in 
energy from |m) by hui nrn and tunnels to the voltage state 
with a rate r n . (c) The photograph shows a Nb/A10 x /Nb 
device. Not seen is an identical SQUID coupled to this device 
intended for two-qubit experiments; the second SQUID was 
kept unbiased throughout the course of this work. 



II. dc SQUID PHASE QUBIT 

Figure [TJa) shows the circuit schematic for our dc 
SQUID phase qubiti^ 3 - The qubit junction Jl (with crit- 
ical current Iqi and capacitance C\) is shown on the 
left. It is isolated from the current bias source lb by 
geometrical inductances L\ and L 2 and the second junc- 
tion J2 (with Iq2 and C 2 )- In order to independently 



control the currents in the two arms of the resulting 
dc SQUID, a current source J/ applies a flux <I> a to the 
SQUID loop through mutual inductance M. Good iso- 
lation of the qubit junction is obtained when L\/M 3> 1 
and Li/ (L-2 + Lj 2 ) » 1, where Lj 2 is the Josephson 
inductance of the isolation junction^ 

For arbitrary values of the bias current lb and flux cur- 
rent If , the dynamics of a dc SQUID can be described by 
2 degrees of freedom corresponding to the phase differ- 
ences across each of the junctions. We, however, operate 
the device by increasing Jf, by Alb while simultaneously 
increasing // by LiAIb/M. This nominally keeps the to- 
tal current through the isolation junction J2 near zero, 
so that the qubit junction current is roughly Further- 
more, one obtains a weak dynamical coupling between 
the junctions by choosing L\ to be large and biasing the 
SQUID so that the two junctions are well out of reso- 
nance with each other^ In this case, the dynamics of 
the phase difference 71 across the qubit junction are gov- 
erned to a good approximation 3 -^ by the Hamiltonian of 
a single current-biased junction ) 35 ' 36 

4E C 2 / h - Irf COS UJ rf t \ 

T~L = —^2~Pi ~ Ej I cos 71 + J — (1) 

Here, Ec = e 2 /2Ci and Ej = I i$ q /2tt are the charging 
energy and Josephson coupling energy of the qubit junc- 
tion, pi = ($o/27r) C±ji is the momentum conjugate to 
71 , and I r f is the amplitude of a microwave drive current 
of frequency w r f. In quantizing 7i, 71 and p\ become 
operators with [71, pi] = i%. 

The second term on the right-hand side of Eq. ((T|) 
defines a one-dimensional tilted washboard potential U, 
sketched in Fig. QJb) . Each well is characterized by the 
classical plasma frequency u p and barrier height At/, 
which can also be expressed as the dimensionless quan- 
tity N s = AU/hojp. A single potential well supports 
roughly A^ s metastable states \n), where the ground state 
|0) and first excited state |1) serve as the basis for quan- 
tum computation^ Motivated by this, the Hamiltonian 
can be expressed in a discrete representation as 

N-l AT-l 
Hn = fr^On W) (n \ + Mlnm \ n ) ( m \ COSLU r ft, 

n—0 n,m— 

(2) 

where TV is the number of states in the well being con- 
sidered, huj nm is the energy-level spacing between states 
\n) and \m) (tunable through lb), and 

nm = — ^ % (n I 71 I rn) (3) 

Z7T n 

is a bare Rabi frequency. 

The simplified Hamiltonian in Eq. ^ neglects tunnel- 
ing through the potential barrier ) 36 ' 37 ' 38 ' 39 a process cor- 
responding to the SQUID spontaneously switching from 
the supercurrent to the voltage state; tunneling from \n) 
occurs at an escape rate T n . As r n+ i/r„ ~ 10 3 for val- 
ues of that we studied, tunneling can be exploited to 
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perform state readout. Note also typically T n <C to n ^ n +i 
for the lowest few levels; for states near the barrier, this 
is not true. The total escape rate of the system to the 
voltage state is 

N-l 

r = V PnnTn, (4) 

Ptot n=0 

where p nn is the occupation probability of state \n) and 
Ptot — X) Pnn is the probability of the system being in the 
supercurrent state; all of the quantities in the equation 
are time dependent during a bias ramp. 

This description of the qubit makes several simplifica- 
tions. For one, we have ignored the quantum states of the 
isolation junction J2.— For typical bias conditions, we 
keep the current through this junction small and its first 
excited state is two to three times higher in energy than 
that of the qubit. However, the higher excited states of 
the qubit junction may come into resonance with states 
of the isolation junction. This approach further ignores 
the role of J2 in determining the bias conditions. Quan- 
tum mechanical simulations of the SQUID show that the 
current through the qubit junction weakly depends on its 
quantum stated shifts in the qubit current are predicted 
to be less than 5 nA (much smaller than Toi ~ 18 pA) for 
our situation, so we have not considered this correction 
below. As a result, the junction parameters in Eq. ([T]) 
need not be equal to the actual values for the qubit junc- 
tion. 

The simulations described in the sections that follow 
require values of T„, oj nmi and (n | 71 | m), all of which 
are specified by lb, I01, and C\. In the absence of dissi- 
pation, analytical expressions for these quantities can be 
obtained by applying perturbation theory to the cubic 
approximation of the tilted washboard potential^ For 
the escape rates and energy levels, we instead use nu- 
merical solutions to the exact potential, which are more 
reliable for states near the top of the barrier. These solu- 
tions come from two methods that give consistent results: 
(i) solving the Schrodinger equation with transmission 
boundary conditions^ and (ii) complex scaling^ 2 - 

We performed experiments on a dc SQUID phase 
qubit fabricated by Hypres, Inc42 using a 30 A/cm 2 
Nb/A10 x /Nb trilayer process [see Fig. [Tfc)]. The qubit 
and isolation junctions had areas of 10 x 10 pm 2 and 
7x7 pm 2 , respectively, and the inductances were roughly 
Li = 3.4 11H, L 2 = 30 pH, and M = 13.4 pH, 44 The error 
in these and other fit parameters reported in this paper 
is about one unit in the least significant digit. We have 
not attempted a thorough analysis of the uncertainties 
in all the parameters due to the complexities of the non- 
linear functions involved. The microwave current I r f was 
carried inside the refrigerator on a single length of lossy 
stainless coax. It was coupled to the qubit via an on- 
chip 1 fF capacitor C r f [see Fig. Q^a)]; while the large 
impedance mismatch this produced potentially improved 
the isolation, it did make an independent calibration of 
the power reaching the qubit difficult. The device was 



mounted in a superconducting aluminum box that was 
attached to the mixing chamber of a dilution refrigera- 
tor with a base temperature of 20 mK. The refrigerator, 
located in an rf shielded room, was surrounded by a mu- 
metal cylinder. In addition, the measurement and bias 
lines were protected from noise at room temperature by 
discrete LC filters and copper powder filters at the mix- 
ing chamber. 

The escape rate of the qubit junction was measured by 
simultaneously increasing and as described above, 
while monitoring the voltage on the current bias line 
of the SQUID. We measured the time interval between 
the start of the ramps and when the system tunneled to 
the voltage state. Repeating this procedure many times 
(~ 10 5 ) at a rate of about 250 Hz yielded a histogram 
of switching times, from which the escape rate could be 
calculated.— Decreasing the repetition rate did not yield 
a significant change in the escape rate, suggesting that 
heating due to the device switching to the voltage state 
had a minimal impact on the measurements. Because the 
large loop inductance of the SQUID resulted in about 20 
possible flux states, the device was initialized to the zero 
trapped flux state with a flux shaking procedure before 
each repetition of the measurement cycled We typically 
applied 50 oscillations of the flux current, which yielded 
a success rate better than 98%. 

Figure Ufa) shows a spectrum of transitions for the 
qubit junction. With a microwave drive of fixed fre- 
quency applied to the device, the total escape rate was 
measured while ramping the biases. The open circles in- 
dicate the values of lb where the microwave drive caused 
the largest enhancement in the escape rate (over its back- 
ground value in the absence of microwaves). 

As shown by the solid line, we fit these points to the 
theoretical values of woi, yielding ioi = 17.930 ^A and 
C\ = 4.50 pF. These values are likely to differ from 
the actual critical current and capacitance of the qubit 
junction, both in principle, due to the simplifications in 
the Hamiltonian mentioned above, and in practice, as 
this fitting procedure is sensitive to inaccuracies in the 
simultaneous ramping of and If, constant offset flux 
that biases the SQUID (from, for example, trapped flux 
near the device), and drifts in the detection electron- 
ics. Nonetheless, when we repeated the measurement 
at elevated temperature, the location of the transitions 
between higher levels agreed to within 0.2% of predic- 
tions for LO12 (dashed) and U23 (dotted) obtained with 
the effective parameters^ 2 - In addition, two-photon tran- 
sitions between these levels occurred at high power at 
the expected frequencies. Multiphoton transitions be- 
tween the ground and first excited states of a variety 
of junction qubits have previously been observed and 
characterizedj 19 ! 45 ' 46 Over the course of three months 
(during which the refrigerator remained near its base 
temperature), the junction parameters found from the 
spectral fits varied by roughly 1%. Thus, we repeated 
the spectroscopy measurement often to obtain the values 
of Iqi and C\ needed for the simulations; these values are 
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17.69 17.73 17.77 17.81 

FIG. 2: Qubit energy-level spectroscopy and tunneling escape 
rates, (a) Open circles show the resonance frequency of the 
transition between the ground and first excited states of the 
qubit, measured at 20 mK. The scatter in the values is indica- 
tive of the uncertainty in the measurement. Also plotted are 
theoretical values of luqi (solid), W12 (dashed), and W23 (dot- 
ted) for J01 = 17.930 fj,A and d = 4.50 pF. (b) At high bias, 
the measured background escape rate (open circles) agrees 
with the predicted ground-state escape rate To (solid line) for 
the junction parameters given above. Calculated Vi, r%, and 
T3 are plotted as dashed, dotted, and dashed-dotted lines. In 
both plots, the bottom axes show the total bias current lb, 
while the top axes indicate the normalized barrier height N s , 
calculated using the extracted junction parameters. 



listed in the caption of each figure. 

The open circles of Fig. [2Kb) show the measured back- 
ground escape rate in the absence of microwaves. The 
solid line indicates the theoretical value of Tq, calculated 
with the values of Iqi and C\ extracted from the spec- 
trum. At low bias current, the measured escape rate ex- 
ceeds Tq, suggesting the presence of excited state popula- 
tion even though the refrigerator was at 20 mK. Based on 
experiments described elsewhere^ we believe that fea- 
tures such as the peak near I b = 17.73 /iA are due to 
population in |2) generated when noise on the leads at 
the resonant frequency of the isolation junction matches 
the 0^2 transition frequency of the qubit. At high bias, 
the escape rates exceed the excitation rates of the noise 
and the total escape rate converges to To; a similar effect 
occurs for thermal excitations^ 7 - While noise complicates 



the situation, the overall behavior suggests that both <x>oi 
and Tq of the dc SQUID are described by the same one- 
dimensional tilted washboard potential. Also shown in 
Fig. [2Jb) are numerical predictions for Ti (dashed), T2 
(dotted), and T 3 (dashed-dotted). 

We observed Rabi oscillations between the ground and 
first excited states by turning the microwave current I r f 
on when If, was at the value where loqi = ^o r f and mea- 
suring the time-resolved escape rate while I r f remained 
at a constant value. This serves as a simple method 
of monitoring the evolution of the state populations. 7 
Although the biases continued to increase during the 
Rabi oscillations, the ramp rates were reduced (with 
dlt/dt ss 0.01 A/s) before I r f was turned on, so that the 
level spacing luqi changed by a negligible amount dur- 
ing the escape rate measurement. The symbols in Fig. [3] 
show r due to a 6.2 GHz microwave drive for a range 
of powers Ps at the microwave generator. As expected, 
the oscillation frequency increases with power and deco- 
herence causes the amplitude of the oscillations to decay 
with time. 

However, there are two unexpected features in the 
data. First, the escape rate at long time (once the os- 
cillations have decayed away) increases with Ps , whereas 
for an ideal two-level system, the excited state population 
saturates at high power. Furthermore, the highest mea- 
sured escape rates far exceed the value of 2.2 x 10 6 1/s 
predicted for Ti at this bias current, strongly suggesting 
that the states |2) and higher are becoming occupied. As 
the escape rates from these levels are very large, only a 
small population would be required to account for the ob- 
served r. Second, at high power, the oscillation maxima 
increase over the first few cycles [see Fig. [3ja)]; this is 
due to the rise time of the microwave current amplitude, 
which will be discussed in Sec. llVl 

In Sec. IIII| we will show how more complete measure- 
ments reveal that the higher excited states impact the 
Rabi oscillation frequencies in ways that are in quantita- 
tive agreement with predictions from the Hamiltonian of 
Eq. but are difficult to see in Fig.[3J As shown by the 
solid lines in the figure, most of the features of the mea- 
sured escape rate are captured by a multilevel density- 
matrix simulation which will be described in Sec. IIVI 



III. DETUNING AND STRONG FIELD 
EFFECTS 

We first model Rabi oscillations using the rotating 
wave approximation, in the absence of tunneling and 
dissipation. This provides a way of predicting the os- 
cillation frequency for a wide range of experimental 
parameters j 9 ' 12 i 27 ' 28 i 29 ' 48 In the rotating frame corre- 
sponding to the drive frequency uj r f luqi, the Hamilto- 
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FIG. 3: Rabi oscillations in the escape rate Y were induced 
at lb — 17.746 fiA by switching on a microwave current at 
t — with a frequency of 6.2 GHz (resonant with the — > 1 
transition) and source powers Ps between —12 and —32 dBm, 
as labeled. The measurements were taken at 20 mK. The 
solid lines are from a five-level density-matrix simulation with 
Joi = 17.930 fiA, Ci = 4.50 pF, T\ = 17 ns, and T$ = 16 ns. 

nian of Eq. ((2|) simplifies tc 12 

N-l -, N-l 

H§ w =J2 nA n\n)(n\ + - fM' nm \n)(m\, (5) 

n— n,m— 

where A n = uj$ n — ncj r f and 

^nm = ^nrn ^m— n+2s— 1 [ J (6) 

for n < to and f2j„ n = Here, J n (x) is the nth order 

Bessel function and VL nm is defined by Eq. 

Differences between the N eigenvalues of this Hamil- 
tonian specify effective multilevel Rabi frequencies or 
modes of the system. We find it convenient to label 
these frequency differences by the states |n) with the 
largest weight in the two corresponding eigenfunctions. 
For the parameter regime of interest here, the differences 
can be uniquely classified by two states, which we denote 
by ^R,nm- For example, il/j,02 denotes the Rabi oscilla- 
tion frequency between eigenstates that describe a two- 
photon transition between |0) and |2). While approxi- 
mate analytical solutions can be obtained for a three-level 
junction system^ we numerically found the eigenvalues 
of the rotating wave Hamiltonian for systems with up to 
seven levels in the simulations that follow. 



FIG. 4: Rabi oscillation frequency Q^oi at fixed bias as a 
function of microwave current I r j. Extracted values from data 
(including the plots in Fig. [3j are shown as circles, while the 
rotating wave solution is shown for two- (dashed line) and five- 
(solid) level simulations, calculated using ioi = 17.930 piA and 
Ci = 4.50 pF with Lu r f/2ir = 6.2 GHz. 

We fit the escape rates in Fig. [3] (and additional data 
for other powers not shown) to a decaying sinusoid with 
an offset. The extracted frequencies are shown with cir- 
cles in Fig. [U To compare to theory, Hr qi, calculated 
using the rotating wave solution for a system with five 
levels, is shown with a solid line. The implied assump- 
tion that the oscillation frequencies of T and p\\ are 
equal, even at high power in a multilevel system, will 
be addressed in Sec. [V] fn plotting the data, we have 
introduced a single fitting parameter 117 nA/VmW that 
converts the power Pg at the microwave source to the 
current amplitude I r f at the qubit. Good agreement is 
found over the full range of power. 

As I r f increases in Fig. |4j the oscillation frequency is 
smaller than the expected linear relationship for a two- 
level system (dashed line). This effect is a hallmark of a 
multilevel system and has been previously observed in 
a similar phase qubit There are two distinct phe- 
nomena that affect — » 1 Rabi oscillations in such a 
device i 12 i 26 ' 28 i 30 To describe these, we must first define 
resonance as occurring when the Rabi frequency rJ^oi is 
at a minimum, as the detuning between the microwave 
drive frequency w r f and level spacing loqi is varied (for a 
fixed drive power). In a two-level system, this happens 
for ojoi = Urf ■ For our phase qubit, resonance occurs 
when u)qi < u r f due to the decreasing level spacings with 
increasing state | n) [see Fig. Ufa)]. This shift in resonance 
is an analog of the ac Stark effect. Resonance shifts occur 
under strong driving in other superconducting qubits as 
well^ 

However, this shift does not explain the data shown in 
Fig.[3] which were taken at a fixed bias current (and thus 
off resonance at high power). The higher levels also affect 
the frequency of the Rabi oscillations. We will refer to 
the minimum value of SIr nm as the on-resonance Rabi 
frequency ri™™ m . The suppression of below S7 i 

leads to the effect seen in Fig. 2J 
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Both of these effects become significant as fioi ap- 
proaches luqi — u>i2, which is a measure of the anhar- 
monicity of the system ; 12 i 28 i 29 in the case of Fig. |4j 
u>oi/27r = 6.2 GHz and lui2/2tt = 5.5 GHz. Clearly, these 
shifts need to be considered when working at high power 
or at low current bias. 

In order to follow experimentally the shift of the res- 
onance condition, it was necessary to measure Rabi os- 
cillations for different detunings of the microwave drive. 
We chose to do this by keeping the drive frequency oj r f 
fixed and changing the level spacing luqi (through if,), be- 
cause the power transmitted by the microwave lines had 
a nontrivial frequency dependence. Figure [5ja) shows 
a grayscale plot of Rabi oscillations measured from such 
an experiment, where black represents a high escape rate. 
Each horizontal line is the escape rate versus time due 
to a microwave current of 6.5 GHz and —11 dBm, which 
was turned on at the value of the current bias lb indicated 
on the vertical axis. While the measurements were per- 
formed at 110 mK, this is not expected to have a signifi- 
cant impact on the Rabi oscillations, as the temperature 
was well below fiwoi/fce ~ 325 mK. 13 

From Fig. EJa), we see that the oscillation frequency 
depends on the current bias. This variation can be seen 
more readily in Fig.[5jb), which shows the power spectral 
density (calculated as the absolute square of the discrete 
Fourier transform) of the escape rate data in Fig. E(a) 
between t = 1 and 45 ns. For this plot, each horizontal 
line has been normalized to its maximum value (black) 
in order to emphasize the location of the dominant fre- 
quency. Three distinct bands are visible. 

For this data set, the level spacing ujoi/2it is equal 
to the microwave frequency uj r f/2n = 6.5 GHz at 
lb = 17.614 /iA. The band with the highest cur- 
rent in Fig. EJb) is centered about = 17.624 /iA, 
suggesting that — > 1 Rabi oscillations are the dom- 
inant process near this bias. For slightly higher or 
lower the oscillation frequency increases as Qr,oi ~ 

\J + (u> r f — woi) 2 , in agreement with simple two-level 
Rabi theory, leading to the curved band in the grayscale 
plot. 

The other bands correspond to oscillations between the 
ground state and higher excited states. Due to the an- 
harmonic level structure, wo2/47r is 6.5 GHz at a smaller 
current bias lb = 17.594 /iA; thus a second band appears 
there, corresponding to two-photon — > 2 Rabi oscilla- 
tions. Similarly, three-photon — > 3 oscillations are vis- 
ible near 17.572 //A where cuq3/6tt = 6.5 GHz. Finally, 
large escape rates occur near 17.549 /iA corresponding to 
a four-photon — > 4 transition, although no oscillations 
are apparent in Fig. [U(a). Note that the effective junc- 
tion parameters /oi and C\ (given in the caption) used to 
predict u> nm and other level properties are slightly differ- 
ent than those for Figs. HH4] as the data sets were taken 
two months apart. 

The rotating wave solution provides a simple way to 
predict the oscillation frequencies. Calculations of f^on 
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FIG. 5: (Color online) Multiphoton, multilevel Rabi oscilla- 
tions plotted in the time and frequency domains, (a) The 
escape rate V (measured at 110 mK) is plotted as a func- 
tion of the time after which a 6.5 GHz, —11 dBm microwave 
drive was turned on and the current bias /j, of the qubit; T 
ranges from (white) to 3 x 10 s 1/s (black), (b) The nor- 
malized power spectral density of the time-domain data from 
t = 1 to 45 ns is shown with a grayscale plot. The dashed 
line segments indicate the Rabi frequencies obtained from the 
rotating wave model for transitions involving (from top to 
bottom) 1, 2, 3, and 4 photons, evaluated with junction pa- 
rameters ioi = 17.828 /iA and C\ — 4.52 pF, and microwave 
current I r f = 24.4 nA. Corresponding grayscale plots calcu- 
lated with a seven-level density-matrix simulation are shown 
in (c) and (d). 



using Eq. ([5]) for a seven-level system are shown as 
dashed curves in Fig. E^b) for the four lowest multipho- 
ton transitions (n — 1, 2, 3, 4). The microwave amplitude 
I r f = 24.4 nA is the only free parameter in the calcula- 
tion and the rotating wave solution using this value repro- 
duces the oscillation frequencies of the different processes 
well, even at large detuning. 

Figure [Hh) shows that the minimum oscillation fre- 
quency n]f 'o\/2ir = 540 MHz of the first (experimental) 
band occurs at lb = 17.624 /iA, for which uiqi/2it = 
6.4 GHz. This again indicates an ac Stark shift of this 
transition, which we denote by Awoi = &rf — c^oi ~ 
2ir x 100 MHz. In addition, the higher levels have sup- 
pressed the oscillation frequency below the bare Rabi fre- 
quency of fi i/27T = 620 MHz [calculated with Eq. ©]. 

We repeated this analysis of resonance for five differ- 
ent microwave powers, with data taken at a later date. 
Figure [6] shows experimental results for the on-resonance 
Rabi frequencies ^i^™ m and Stark shifts Auj nm with the 
corresponding results from a five-level rotating wave so- 
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FIG. 6: (Color online) The (a) on-resonance Rabi oscilla- 
tion frequencies fi^'oi and f2™02 and (b) resonance frequency 
shifts AtJoi = Urf — Woi and Ai^o2 = 2uj r f — luq2 are plotted as 
a function of the microwave current, for data taken at 110 mK 
with a microwave drive of frequency lo t j /2n = 6.5 GHz and 
powers Ps — —23, —20, —17, —15, —10 dBm. Values extracted 
from data for the — > 1 (0 — > 2) transition are plotted as open 
circles (filled squares), while five-level rotating wave solutions 
for a junction with Ioi = 17.736 fiA and Ci = 4.49 pF are 
shown as solid (dashed) lines. In (a), the dotted line is from 
a simulation of a two-level system. 

lution for the — > 1 (circles for data and solid lines for 
theory) and two-photon 0^2 (squares and dashed lines) 
transitions. Here, the power calibration is 84 nA/VmW 
for oj r f/2n = 6.5 GHz. Figure \@[ a) differs from Fig. 2J 
because in the former, if, was varied at each power to 
give the minimum oscillation frequency; by staying on 
resonance in this way, the effect of the higher levels on 
the — ► 1 oscillations is maximized. 

The resonant oscillation frequencies fij^ o\ and fijjj 02 
in Fig. [SJ a) are well described by Eq. ([5]) over the full 
range of I r f . The deviation between the — > 1 oscillation 
frequency and the values expected in a two-level system 
(dotted line) increases with I r f. Similar measurements 
taken at a higher If, show a smaller frequency suppression 
over a similar range of I r fr^ as expected for a system 
with stronger anharmonicity. 

The resonance shifts Ac^oi and Awo2 m Fig.^b), how- 
ever, differ significantly from the model predictions. As 
a change of 10 MHz in ujoi/2it corresponds to a roughly 
1 nA change in If,, the discrepancy is difficult to see in 
Fig. [S](b) . The indicated uncertainty in the experimen- 



tal points in Fig. [6]Jb) is roughly 5 MHz due to errors in 
the calibration of cjoi {h)- In addition the uncertainty 
is somewhat larger at low power, where the relatively 
small total escape rates result in poor counting statistics, 
and high power, where the weak dependence on detuning 
makes it difficult to identify the resonant level spacing. 
While the general trend of the shifts is consistent with the 
model, the scatter in the data is large and further work 
would be needed to determine if there are true deviations 
from the multilevel theory. 



IV. DECOHERENCE 

We now model the time dependence of the escape rate 
measured in the experiment. In order to do this, several 
additions have to be made to the treatment of the system 
given in Sec. Hilt the effects of tunneling, other sources 
of damping and noise, and experimental limitations. The 
density-matrix formalis m 51 ' 52 provides a straightforward 
scheme for including nonunitary processes and has been 
previously applied to the lowest two or three levels of the 
phase qubiti£ ' 28 i 29 ' 30 i 48 ' 53 For a system with N levels, we 
assume that the evolution of the qubit's reduced density 
matrix p is given by the modified Liouville-von Neumann 
equation 

^ = ~[H N ,p]-Gp-Rp-Dp. (7) 

Here, we use the discrete Hamiltonian in Eq. |(3J), so 
that the diagonal elements p nn give the occupancy of 
the states \n). Tunneling is characterized by the tensor 
G, where [Gp] nm = (T n + T m ) p nm /2, which leads to a 
decay of all elements of p. The tensors R and D account 
for two distinct decoherence mechanisms. 

The Bloch-Redfield tensor R models the effects of the 
system being in equilibrium with a thermal bath at tem- 
perature TM- This tensor leads to the decay of the diago- 
nal elements of p (dissipation) as well as the off-diagonal 
elements (decoherence). The coupling to the bath, as- 
sumed linear in 71, is parametrized by R\ (to), which is 
the inverse of the real part of the total admittance that 
shunts the qubit junction, evaluated at angular frequency 
w. In Eq. 0, 

JV-l 

[ R P\nm = X! RnmklPkl, (8) 
k,l=Q 

where 

AT-l JV-1 

Rri7nkl — 'Ylmnk "iknml ~t~ $lm ^ ^ Tnrrk "F ^nk ^ ^ ^/mrrl 

r=0 r=0 

(9) 

and 
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(I | 71 1 m) (n | 71 1 fc) 



(1 — <5nfc) W Tl fe- 



exp i 



hu n kSgn(n — k) 
2k„T 



sinh (hw nk /2k B T) 



$nk 



2k B T 



(10) 



In the following, we assume that R\ is independent of 
frequency, in which case we can define a dissipation time 
Ti = R\C\. In this limit, Eq. ([5]) gives well known inter- 
level transition rates i 37 ' 38 i 55 For example, thermal exci- 
tation from \n) to |m) (with m > n) occurs at a rate 



(a) 



w + = 



27n mm 

2 /£o 
ft V 2tt 



|(ra|7i \m)\* 



Ri J exp (hw nm /k B T) - 1 



(.11) 



As required by detailed balance, decay from \m) to \n) 
occurs at a rate W~ m = W+ n ex P (f^Jnm/k B T). In addi- 
tion, Eq. ([5]) specifies the decoherence due to this dissipa- 
tion. Given uo nm and the matrix elements of 71 , all of the 
thermal rates are specified by T\, which is set to 17 ns for 
the simulations below. This value comes from additional 
measurements of the escape rate of the device in the ab- 
sence of microwaves over a range of temperatures. 4 ? 

We find that this treatment of dissipation alone is 
insufficient to capture the decay of Rabi oscillations, 
suggesting that an additional decoherence mechanism is 
present. This we model by the tensor D in Eq. ([7]), which 
has the form 



D P = Xn ( L "P L n ~ ~ L tL, 



P - 



T>P L \ L i 



(12) 



where L n are Lindblad operators with strengths A„ . The 
best overall agreement with the measurements (see fur- 
ther discussion below) is found for a set of operators 
L n = \n) (n\, where A n = l/T^ and n ranges from 
to N — 1. This leads to an exponential decay of each 
of the off-diagonal elements of p with a common time 
constant given by the dephasing time and no change 
in the diagonal elements. In the simulations of this 
section, we set = 16 ns, yielding a coherence time 

T 2 = [1/ (271) + l/T^P 1 - 10-9 ns. 

For N = 2, Eq. (0 reduces to the optical Bloch 
equations ; 51 ^ 2 for which analytical solutions exist. How- 
ever, for N > 2 the master equation is easily numerically 
integrated to obtain the time dependence of p, without 
making the rotating wave approximation. With the state 
occupation probabilities p nn in hand, the total escape 
rate can be calculated with Eq. (|4]) and compared with 
experiment. 

As mentioned earlier, one of the more striking features 
of the series of — > 1 Rabi oscillations in Fig. [3] is that 
Too, which we define as the steady escape rate the system 
approaches as the oscillations decay away, increases over 
the full range of measured power. The circles in Fig. [7] 
show experimental values of as a function of the bare 
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FIG. 7: Long-time escape rate Too as a function of the 
bare Rabi frequency fioi- ( a ) Too is plotted for a density- 
matrix simulation of a system with two (dashed-dotted) , 
three (dashed), four (dotted), and five (solid) levels, I01 = 
17.930 pA, and C*i = 4.50 pF. (b) at low power (for a five- 
level system) is affected by the inclusion of a small microwave 
current at W02 (dashed) and inhomogeneous broadening in ad- 
dition to this noise current (solid). One set of experimental 
data, for source powers between —68 and —11 dBm at 6.2 
GHz, is plotted in both panels with open circles. 



Rabi frequency f^oij which was calculated from the mi- 
crowave source power by using the fit in Fig. [3J We will 
now use the density-matrix simulations to understand 
these escape rates. 

The simulations suggest that T^, is not sensitive to 
the time evolution of p, and is only weakly dependent 
on T\ and at high power. However, Y ^ does depend 
on the number of states that participate in the dynamics 
and their individual escape rates T n . It should be noted 
that the actual populations do not reach steady state, 
but continuously decay due to tunneling. The dashed- 
dotted line in Fig. [Tfa) shows the calculated for a 
two-level system (N — 2). It increases up to a value of 
Ti/2 at roughly Sloi = l/V^T^T; the data show only a 
subtle shoulder, masked by an overall steady increase, 
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near this escape rate. Results for three (dashed), four 
(dotted), and five (solid) levels are also plotted in the 
figure, displaying increasing agreement with experiment. 

This agreement between theory and experiment at high 
power is an indication that transitions to the higher 
states, or leakage in the context of quantum computation, 
are occurring as expected. As power increases, the occu- 
pation of higher excited states increases and this effect is 
magnified by their larger escape rates. For example, with 
the five- level simulation at f2oi/27r = 1 GHz the state |3) 
has a 2% occupation probability, but accounts for 60% of 
Too. Note that microwave pulse shaping will not reduce 
this effect. That is, while leakage can be minimized at the 
end of a single qubit operation, higher states are always 
populated during the pulse i 27 ' 28 > 29 Such pulses will have 
enhanced escape rates comparable to those seen here for 
a constant I r f. 

These simulations do not do a good job of explain- 
ing the value of at very low power. The reason for 
this can be seen in Fig. [2fb), where the measured es- 
cape rate in the absence of microwaves exceeds the pre- 
dicted value of Tq. In particular, at the bias current 
of If, = 17.746 /iA where the Rabi oscillations were per- 
formed, the measured escape rate is ~ 2 x 10 4 1/s. While 
this excess escape rate could be attributed to a thermal 
population at T = 56 mK (whereas the refrigerator ther- 
mometer read 20 mK), raising the temperature of the 
simulation leads to a significant enhancement in Too that 
extends up to moderate microwave power, an effect not 
seen in the data. Instead, separate measurements^ 4 - in- 
dicate that the spurious features in the measured T are 
largely due to a population in |2), beyond that expected 
from a thermal bath at 20 mK. Assuming that this was 
a result of noise on the bias lines, we included in the sim- 
ulation an additional microwave source at a frequency of 
ujq2', a microwave current amplitude of 0.2 nA reproduces 
the background escape rate measured at lb — 17.746 /iA 
(corresponding to Figs. El HI and[7|) and leads to a popula- 
tion P22 = 3 x 10 -5 . The coherent effects of this drive are 
negligible, as Jlo2/27r = 0.5 MHz, which is much smaller 
than 1/Ti or l/T^. The results of a five-level simulation 
with the extra source included are shown with a dashed 
line in Fig.[7jb). Compared to the simpler simulations in 
Fig. Elja), they show improved agreement at the lowest 
powers, with little change above £!oi/27r = 10 MHz. 

To this point, we have not included effects from inho- 
mogeneous broadening, although our spectroscopic mea- 
surements suggest its presence. For = 10.9 ns, which 
describes many of the experiments below, the expected 
full width at half maximum of a resonance peak is roughly 
30 MHz. However, we commonly find peak widths of 50 
MHz at bias currents where tunneling makes a negligible 
contribution. The remaining broadening may be due to 
current noise at frequencies much lower than 1/Ti, whose 
contribution to the width scales with the spectral slope 
diOQi/dl},^ This can be modeled by taking into account 
the frequency content of the noise and using the stochas- 
tic Bloch equations^ Instead, we mimicked the inferred 



spectroscopic broadening simply by running the simula- 
tion for a range of bias currents and then convolving the 
resulting escape rate (at a given time) and a Gaussian 
with standard deviation 07 = 1.5 nA (corresponding to 
a 35 MHz spread in cJoi/27r). 

Calculations of for a five-level system with noise 
at W02 and this inhomogeneous broadening are drawn 
with a solid line in Fig. [7jb). The extra broadening has 
negligible effect at high power (where Rabi oscillations 
are observed), but does bring the simulation into better 
agreement with data near Q,qi/2it =10 MHz. The small 
remaining discrepancy could be due to detuning arising 
from a misidentification of the mean value of If,. For ex- 
ample, performing the simulation at 17.745 rather than 
17.746 /iA results in an underestimate of the measured 
values. Nevertheless, the overall agreement over three 
decades in the Rabi frequency is very good. While the 
calibration of the microwave current I r f came from a fit 
to data, one conversion factor reproduces both the os- 
cillation frequencies (Fig. [4]) and long-time escape rates 

(Fig. ®. 

We next consider the time dependence of the escape 
rate for the data plotted in Fig. [H Here, a 6.2 GHz 
microwave pulse nominally 30 ns long was applied on 
resonance with the — >• 1 transition of the qubit junc- 
tion. The measured escape rate shows Rabi oscillations 
followed by a decay back to the ground state once the 
microwave drive has turned off. This decay appears to 
be governed by three time constants. Nontrivial decays 
have previously been reported in phase qubits 1 ^ and we 
have found them in several of our devices. 

Accurately simulating this experiment requires knowl- 
edge of the time dependence of the microwave pulse, 
which was created by the internal gate of a Hewlett- 
Packard 83731B synthesized source^ without any fur- 
ther filtering. We measured the pulse envelope at the 
source's output using a digital sampling oscilloscope [see 
inset of Fig.[UJa)]. Ignoring any distortion of the pulse be- 
fore it reached the qubit junction, this was taken as the 
microwave amplitude I r f (t) . For the microwave power 
used in Fig. [8l the maximum value of I r f was taken to 
be 11.7 nA (a value obtained from the power calibra- 
tion in Fig. The solid line in Fig. Ufa) shows the 
calculated escape rate for a five-level simulation in the 
presence of noise at UJ02 and inhomogeneous broadening. 
For this time-domain plot and others discussed below, 
we have convolved the simulation and a Gaussian with 
a full width of 150 ps to remove very fast, small oscilla- 
tions (due to highly detuned multiphoton processes) that 
could not be seen in the experiment due to insufficient 
time resolution. 

The main part of the Rabi oscillation is reproduced 
well. In particular, the second oscillation maximum has 
a larger escape rate than the first, due to the 7 ns it 
takes for I r f to reach its maximum. The first part of 
the decay is also reproduced, which from the simulation 
should correspond to the emptying of state |2). However, 
the data also show a slow decay constant longer than 
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FIG. 8: (Color online) Response to a microwave pulse, (a) 
The measured escape rate T (circles), due to a 6.2 GHz, 
—20 dBm microwave pulse nominally 30 ns long, shows Rabi 
oscillations followed by a decay governed by multiple time 
constants. The solid line is the result of a five-level density- 
matrix simulation for a junction with Joi = 17.730 fj,A and 
Ci = 4.46 pF, using the measured microwave pulse amplitude 
(with a maximum of 11.7 nA) shown in the inset. The sim- 
ulation also produces (b) the contribution to the escape rate 
and (c) the normalized occupation for each of the five levels, 
as labeled. 



50 ns that the simulation does not explain. The longer 
time is inconsistent with our thermal measurement of T\ 
and instead may be indicative of the qubit interacting 
with an additional quantum system. It appears, though, 
that this extra degree of freedom does not significantly 
affect the description of the Rabi oscillations. Note that if 
this longer time constant were the dominant relaxation 
process at high power, two-level saturation would have 
occurred at a lower power in Fig. [7fa). 

The density-matrix model can also be used to predict 
the escape rates for the power series in Fig. [3] calculations 
are plotted with solid lines in that figure. The presence 
of the noise signal at W02 has little effect and inhomoge- 
neous broadening decreases the escape rate at the lowest 
microwave power slightly. The maximum value of I r f (t) 
was again calculated for each power using the fit in Fig. 2] 



The T,f, value of 16 ns used in the simulations was chosen 
to best reproduce the decay envelopes over the full range 
of powers. At the highest power, there is a discrepancy 
in the oscillation maxima and minima, perhaps indica- 
tive of inaccurate parameters for levels |3) and |4). The 
envelope of the microwave turn-on is also reflected in the 
shape of the escape rate. 

As a final test of the model, we examine multiphoton 
transitions in the system. Figure H^c) shows a grayscale 
plot of the escape rate calculated with a seven-level 
density-matrix simulation. Nearly all of the features 
seen in the data of Fig. [SJa) are present in the simu- 
lation. As the gray scales are identical, a small quanti- 
tative disagreement is visible, particularly for the 0^3 
three-photon transition. Figure Efd) shows the normal- 
ized power spectral density calculated from Fig. [5](c) ; it 
agrees well with the data in Fig. H^b). 

Figure [5] shows line cuts of the time-domain data (cir- 
cles) and simulation (solid lines) from Fig. [5] at bias cur- 
rents of (a) 17.623, (b) 17.596, and (c) 17.571 fiA. These 
values of lb correspond to the resonances of the one, two, 
and three-photon transitions at the high power at which 
the data were taken. While the decay time and the long- 
time escape rate are reproduced well for the three 
transitions, the first few nanoseconds are not captured 
fully. This could be due to distortion of the microwave 
current by the coupling capacitor C r f or other line mis- 
matches. 

As with all of the simulations discussed in this sec- 
tion so far, the solid lines in Fig. [5] were calculated with 
Ti = 17 ns and = 16 ns. This gives T 2 = 10.9 ns and a 
Rabi decay time roughly equal to the two-level valu o 51 i 57 
of T' = [1/ (2Ti) + 1/ {2T 2 )T 1 = 13.3 ns, consistent with 
measurements of the decay envelope of the escape rate. 
If pure dephasing were not present then T' would be 
22.7 ns, which is significantly longer than what is ob- 
served. Although unlikely, our thermal measurement of 
T\ = 17 ns could be incorrect. To examine this possibil- 
ity, simulations with T\ = 10 ns and Ti — 2T± (which also 
give T 1 = 13.3 ns) are shown with dotted lines in Fig. [9] 
They are nearly identical to the curves calculated with 
dephasing, although is somewhat smaller. Thus this 
data set alone cannot rule out dissipation-limited deco- 
herence. However, the shorter T\ reduces the prediction 
for Too by roughly 15% over the full range of measured 
powers in Fig. suggesting that additional dephasing is 
instead affecting the Rabi oscillations. 

While dephasing is needed to faithfully reproduce fea- 
tures of the experimental measurements, its origin is 
unclear. In the simulations, we assume that each off- 
diagonal term of the density matrix decays with a com- 
mon dephasing time of = 16 ns (along with deco- 
herence due to dissipation, which is transition depen- 
dent). If this dephasing at high power were due to low 
frequency noise (similar to the effects of inhomogeneous 
broadening discussed before) , 55 i 58 the corresponding de- 
phasing operator D in Eq. (fj"2")) could be expressed with 
L = ^On/woi \n) (n| and A = 2/T^ (similar to the 
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FIG. 9: (Color online) Comparison of possible damping sce- 
narios for multi-photon transitions. The measured escape 
rates from Fig. [5Ja) for the (a) — > 1, (b) — > 2, and (c) 
0^3 transitions are plotted with circles. The solid lines 
were calculated with the same simulation (with a common 
dephasing time T^) whose results are shown in Fig.^c). The 
dotted lines were calculated for Ti — 10 ns and no additional 
dephasing, while the dashed lines correspond to Ti = 17 ns 
and pure dephasing with = 16 ns due to low frequency 
noise (see text). 

harmonic oscillator number operator). With this choice 
of operator L and strength A, the — » 1 dephasing is 
unchanged. Simulations with this sort of damping are 
shown with dashed lines in Fig.[5J While the — > 1 Rabi 
oscillations are modeled well, there is far too much deco- 
herence for the higher order transitions. This should also 
be true for other implementations of low frequency noise 
that have dephasing rates that scale as (du> nrn /dlb) 2 ■ 
Thus having the ability to measure a wide range of tran- 
sitions can reveal additional information about decoher- 
ence. 



V. DISCUSSION AND SUMMARY 

We have presented measurements taken on a dc SQUID 
operated in such a way that one of the junctions behaves 
much like a simple current-biased junction. We find that 



the simplified Hamiltonian of Eq. ([T]) gives an accurate 
description of the qubit dynamics. Our observation of 
subtle features of this model provides further confidence 
that it can be applied to design the gates needed for 
quantum computation. The analysis presented here can 
also be extended to describe the behavior of the other 
types of superconducting qubits mentioned in Sec. UJ 

We performed several checks on the model, including 
measuring the Rabi oscillation frequency for one and two- 
photon transitions and comparing the data with predic- 
tions based on the rotating wave Hamiltonian of Eq. 
We find good agreement, despite the model not contain- 
ing any damping. Thus the resonance shifts in multi- 
level Rabi oscillation frequencies are mainly determined 
by the anharmonic level structure and not the T± and Ti 
of the qubit. These multilevel effects were identified un- 
der strong driving where it is clear that they could lead to 
errors in simple two-state rotations. Future experiments 
could mitigate these effects by proper pulse shaping and 
operating at low microwave power. 27 i 29 i 32 

We also used a density-matrix simulation, which in- 
cluded the effects of decoherence, to calculate the tun- 
neling escape rate. The results of the simulation agreed 
well with experimental Rabi oscillation data which were 
acquired with the device biased so that even the ground- 
state escape rate To was measurable. While this simple 
measurement (which ends with the qubit in the finite 
voltage state on every repetition) is not an ideal projec- 
tive measurement to the qubit basis states |0) and |1), it 
is particularly well-suited for the purposes of the current 
work. For one, it is extremely sensitive to excited state 
population due to a microwave drive or thermal transi- 
tions. The leakage information contained in Fig. [7] would 
be considerably more difficult to obtain from the pulsed 
single-shot measurements 9 -^ that excel at determining 
the total population not in the ground state; nonethe- 
less, these techniques have been successfully used to mea- 
sure the second excited state population preciselyj 31 ' 32 i 44 
With the simple tunneling measurement, the large ratio 
r n+ i/r„ provides a natural way of distinguishing which 
excited state is populated. In addition, as nothing is done 
to the qubit to initiate a measurement, this method is 
not limited by a measurement fidelity and does not suffer 
from the potential problems associated with changing the 
bias location in order to perform state readout i 31 i 44 ' 60 i 61 

The main drawback to the escape rate measurement is 
that it does not directly produce individual state popu- 
lations. However, the density-matrix simulations do pro- 
vide each level's contribution to the total escape rate and 
its occupation probability. An example for a — > 1 os- 
cillation is shown in Figs. [UJb) and[8jc). The plots sug- 
gest that for high power, the escape rate during a Rabi 
oscillation is dominated by the contributions p22^2 and 
P33T 3 from states |2) and |3). The figure also shows that 
P22 (t) and P33 (t) have nearly the same form as p\\ (t), 
albeit with a much smaller oscillation amplitude^ Thus 
changes in T are reflective of the underlying oscillation 
of |1), which is why we believe the frequency analysis in 
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Fig. |4] is valid. 

The populations in Fig. [5Jc) have been normalized at 
each time t. While the probability of the junction leaving 
the supercurrent state during a microwave pulse depends 
strongly on I r f and the bias conditions, ptot = 0.42 at 
t = 30 ns for this data set. Thus, while we wish to use 
tunneling as a probe of p nn , it is clearly playing a large 
role in the evolution of the system. However, repeating 
the density-matrix simulation in the absence of tunneling 
(i.e., Gp = 0) results in a less than 5% change in the nor- 
malized populations of |0) and |1) during the oscillation. 
The subsequent decay is significantly different without 
the fast decays due to escape. These insights from the 
simulation hinge on the assumption that tunneling does 
not affect T\ and 7$. 

A few remarks should be made about the many in- 
put parameters required for the simulations. For exam- 
ple, the seven-level density-matrix calculation shown in 
Fig-E^c) needed seven escape rates, six energy- level spac- 
ings, and 28 matrix elements of 71, all of which are func- 
tions of If,. We determined these from the properties 
of a tilted washboard potential defined by 7oi and Ci, 
which were measured independently by low power spec- 
troscopy. To accurately reproduce the measured escape 
rate (particularly at low power), we also had to add sev- 
eral other features to the simulation, each motivated by 
separate measurements. In particular, we added a mi- 
crowave source at W02 to mimic noise (whose magnitude 
was found from the microwave-free escape rate), inho- 
mogeneous broadening due to bias noise (estimated from 
spectroscopic resonance widths), and a finite time resolu- 
tion (consistent with the bandwidth of the detection elec- 
tronics). The dissipation time T\ = 17 ns was estimated 
by an independent thermal escape rate experiment^ 
while the dephasing time 7^ = 16 ns (and the way it was 
incorporated into the simulation) was chosen to maxi- 
mize the agreement with the Rabi data. For simplicity, 
both time constants were assumed to be independent of 
frequency. Removing this condition in the analysis could 
yield additional information about the decoherence in the 
system, as the oscillations are most sensitive to noise at 
the Rabi frequency 58 ' 62 Finally, the conversion between 
the microwave power at its generator and the current 
through the qubit junction was itself calibrated by the ob- 
served oscillation frequency, and the resulting factor was 



allowed to be a function of microwave frequency. Thus a 
wide range of multilevel phenomena was explained with 
just a few truly free parameters. 

In addition, we found that the output of the simula- 
tions (whether that be an oscillation frequency or escape 
rate) converged as the number of levels N was increased. 
For any of the five or seven-level simulations discussed 
here, an additional level produced a change too small to 
be detected by the experiment; more levels were required 
at high power or low lb for the solution to converge sat- 
isfactorily. Surprisingly, the highest included levels are 
predicted to lie above the barrier. However, small errors 
in the energy levels or matrix elements for the highest 
levels do not produce large changes in the final escape 
rates. 

Certain features of the data were not reproduced by 
the simulations. Although not discussed here, the tran- 
sition spectra of our Nb qubits show a series of splittings 
(all less than 10 MHz wide), similar to but smaller than 
those reported in other superconducting qubits ^ 31 ' 63 ' 64 
The extra degrees of freedom responsible for these fea- 
tures were not included in the device Hamiltonian, ex- 
cept perhaps in some effective way through T\ and 7^ 
of the density-matrix simulation. For weak coupling and 
high microwave power, the impact of individual two-level 
systems on Rabi oscillations is expected to be smalL 28 ' 65 
However, it is unclear if a bath of quantum systems could 
be responsible for the anomalous decay in Fig. [5J As 
other groups have observed strong dissipation in Nb tri- 
layer junctions ) 10 ' 13 it is likely that the wiring insulation 
or other details of the fabrication are playing a critical 
role** 6 - While further work is needed to identify (and po- 
tentially eliminate) the actual microscopic sources of de- 
coherence in this system, the multilevel features in the 
dynamics of our qubit are largely understood. 
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